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Abstract. The LISA mission is a space interferometer aiming at the detection of 
gravitational waves in the [10 _4 ,10 _1 ] Hz frequency band. In order to reach the 
gravitational wave detection level, a Time Delay Interferometry (TDI) method must 
be applied to get rid of (most of) the laser frequency noise and optical bench noise. 
This TDI analysis is carried out in terms of the coordinate time corresponding to the 
Barycentric Coordinate Reference System (BCRS), TCB, whereas the data at each of 
the three LISA stations is recorded in terms of each station proper time. We provide 
here the required proper time versus BCRS time transformation. 
We show that the difference in rate of station proper time versus TCB is of the order of 
5 ■ 1CP 8 . The difference between station proper times and TCB exhibits an oscillatory 
trend with a maximum amplitude of about 10~ 3 s. 
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1. Introduction 

The LISA mission [7j is a space interferometer aiming at the detection of 
gravitational waves (GW) in the [~ 10 -4 , ~ 10 _1 ] Hz frequency band. Gravitational 
waves crossing the LISA quasi equilateral triangular constellation are detected through 
the induced change in the station inter-distances. The latter also depend on time due 
to the gravitational field of the Sun [2], mostly, and planets; what we call a "geometry 
(G) effect". 

In order to reach the gravitational wave detection level, a Time Delay Interferometry 
(TDI) method (see [15] for a review and references therein) must be applied to get rid 
of (most of) the laser frequency (LF) noise and optical bench (OB) noise. In other 
words, TDI is needed to bring those LF plus OB (physically indistinguishable) noises 
down to the level of the other noises : quantum (Q), fiber (F), residual proof mass 
motion (PM) noises... The TDI method consists in combining numerically data fluxes 
at the stations (rather than combining the laser beams physically) with an appropriate 
delay. Hence, the so-called TDI observables are closed loop combinations of the different 
laser links with appropriate delays (combination of photon-flight time tij between two 
stations i — > j which correspond to station inter-distances) that cancel (almost all) the 
laser frequency noise and optical bench noise. 

The TDI analysis is carried out in terms of the coordinate time t corresponding to the 
Barycentric Coordinate Reference System (BCRS), the so-called TCB, whereas the data 
at each of the three LISA stations, k — 1,2, 3, is archived in terms of the station proper 
time Tfc. 

We modeled the orbitography of the three stations {xl, Ut = d~x\ /dt with t =TCB and 
i = 1,2,3) classically [6] in the gravitational field of the sole Sun, and the laser links, 
that is the photon-flight time t^, relativistically [2j as a function of the position and 
velocities of stations i and j at emission time. 

We here provide the corresponding analytical proper time versus TCB transformation 
required to apply the TDI analysis. We show that the difference between station proper 
time and TCB reaches about 0.5 s over a one-year mission and exhibits an oscillatory 
trend with a maximal amplitude of 0.0014 s. The difference in rate of station proper 
time versus TCB is of the order of 1.5 ■ 10 _8 s. 

We then show how this proper time versus TCB transformation fits in a general 
relativistic TDI analysis. 

2. Proper time versus coordinate time transformation 

2.1. Classical orbitography 

We assume a classical orbitography for the three LISA spacecraft k — 1,2, 3, in 
the BCRS, as given in [6]. Those BCRS coordinates (x^, yk, Zfc), for arbitrary initial 
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conditions, can be rewritten in terms of rotated keplerian ellipses (x e u y e u k, z e u k) as 
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and uj are the common semi-major axis, eccentricity, inclination and argument of the 
periaster of the three spacecraft orbits, respectively. The optimal inclination of the LISA 
triangle on the ecliptic is v = f + |^ [10] with L = 5 • 10 9 m, the average interferometric 
arm-length. The longitude of the ascending node, Qk, is particular to a given spacecraft 
k and is given in terms of that of the first one with a phase shift 9 k- 

Qk = ^i — &k with 9k = —2 (k 



; 3 



The time parametrization of the orbits is given by the equation of the eccentric anomaly 
^k °f eac h spacecraft, 



^ fc - esin^k 
with the mean anomaly 



(2) 



— (t - t ) + M k0 



in terms of the orbital period, T (provided by Kepler's third law, 2n/T = ^jGM/a 3 ) 
and the mean anomaly of spacecraft k at initial time t , that is Mko = Mk(t = t ). 
Mean anomalies are related to that of the first spacecraft through the phase shift: 



M k = M 1 + e. 



k ■ 



The common spacecraft orbit eccentricity being small, the eccentric anomaly 
equation ([2]) can be developed at 1st order (with respect to e): 

2tt , v „ /2tt 



27T /27T \ 

+ — (t - t ) + Mio + ^ + esin ( — (t - t ) + M 10 + fl fe J . 



(3) 
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BCRS position and eccentric anomaly equations used in [2j correspond to particular 
initial conditions (to = 0, u = 3tt/2, Qi = 3tt/2, M 10 = 0) without any planets (which 
means that both to and Mio are completely arbitrary). When planets are modeled, 
the LISA guiding center, that is the projection of LISA's center of mass on the ecliptic 
plane, is supposed to be about 20 degrees behind the Earth at initial time t . This 
means that t and M 10 are no more arbitrary. For example, if the initial mission time is 
t = on epoch Jan 01 2012 at 00:00:00, we find M w ~ -.955002 from keplerian Earth 
ephemeris for uo = 3it/2, Qi = 37r/2. 

2.2. Relativistic time transformation 

k 

We now wish to compute the time transformation between proper time, r, of the 
clock on board spacecraft k, as a function of coordinate (BCRS) time t, called TCB; since 
coordinate time is the common "language" between the different spacecrafts (k = 1, 2, 3) 
and the time used in the TDI method. 

If we consider only the gravitational field due to the Sun, the time transformation 
is given by (Figure [2]) 

k 2 

ds 2 = c 2 dr : 

\ c z 

-^(i-a-i)* ,4, 

where c is the speed of light in vacuum, Wk = with G, Newton's gravitational 
constant, M, the mass of the Sun, the radial distance relative to the Sun at time t 
and Vk the velocity of spacecraft k at time t in the Barycentric Coordinate Reference 
System. 

We can compute the norm of the satellite keplerian velocity and keplerian radial distance 
using equations ([I]), leading to 

2tt\ 2 j^l + ecos 1 !';; 





r k =a (1 - ecos^ fc ) . (6) 
Substituting those expressions in the time transformation (T4j), we obtain 

c 2 a i l-ecos* fc 2c 2 V T / J 1 — e cos ^ k 



GM T a 2 2tt t a 2 2vr 

~ ; ^ k n — ^k n — esm^k + cst 

c 2 a 2vr * 2c 2 T 2c 2 T 

k , V GMa . . 

~ T ° ~ 2&~ ^ ^ k ~ k0 ^ + 6 k ~ Sm k0 ^ ' 



(7) 



k k 

The integration constant is given by r=r and = ^k[to) = ^fco at initial 
coordinate time t = to- The integration was performed through the change in variable 
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dt = T (1 — ecos 1 !^) /27r <i\l/fc, using the time derivative of the exact expression for the 
eccentric anomaly (j2J), and Kepler's third law. 

From the implicit equation (J2j) providing \l/fc(£), one can compute A as a function of 
time t (Figures [3] and Hj). Alternatively, one could use the 1st order expression (j3J. 

To obtain the figures in section [6J the initial time was chosen as to = 0, the initial 

12 3 

offsets of the clocks where r = 0.1 s, r = 0.2 s, r = 0.3 s with initial conditions 

k 

to — 0, co> = 37r/2, Q 1 = 3n/2 and M w = 0. The cumulative r — t effect reaches about 
half a second over one year. 

In Figure HJ the linear trend is removed. The amplitude of the A oscillatory behavior 
is ~ 10~ 3 s. 

3. LISA model and TDI technique: a coherent general relativistic approach 

k 

We now recall how this proper time of spacecraft k versus TCB transformation, A, fits 
in the TDI picture. Let us consider data flow variables ij, collected in a vector, 

e = (^31) E 12i ^23! ^21; £32> ^13 i e 22i £33) • (8) 

We briefly describe the content of the above vectojjjl using Figure [5j 
Each fiber data-flow variable, 

£jj (t) = (T^(t)-r k ^t))/2 with j = 1,2, 3, (9) 

consists in an antisymmetric combination of the two data flows at station j, 

T ijk (t) = T ijk (10) 

w t, at j with LF,OB,F,PM noises v ' 

and T k ^(t), that compare the frequencies of the two local lasers (jk and ji) via fiber 
links. 

Each laser-link data-flow variable, 

eij(t) = EijA ki {t)-A ij {t) with i ^ j,j k,k i and i, j, k = 1,2, 3, (11) 

consists in comparing, at station j, the frequency shift between the incoming laser beam 
i — >• j and the local laser aligned in direction j — > i at station j, 

A l >(t) = A 11 , (12) 

v ' t, at j with LF,Q,OB,PM noises + GW signal + G v ' 

with that at station i, A kl (t), considering that the light has travelled a path i — > j which 
translates into a corresponding time delay tij. In terms of operators, it means dealing 
with the time-delay operator so that for a given function of time / (t) , 

^(/(t))=/( t _ty) . (13) 

The (coordinate-)time transfer between spacecraft i and spacecraft j, 

tij = t rec j t em j , (14) 

§ Our notations (£31, £12, £23; £21, £32, £13; £11, £22, £33) correspond to (Ui, U 2 , U 3 ; -V\, -V 2 , -V3', 
Zi,Z 2 ,Z 3 ) in reference [8j 
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is the barycentric coordinate night-time of photons between spacecraft % (photon emitted 
at time t em j by spacecraft i) and spacecraft j (photon received at time t rec j by spacecraft 
j). This quantity is useful in TDI techniques: data-flow variables are combined with 
appropriate time delays to form TDI observables which are (almost) free of LF and OB 
noises: 



rn 



in a x 



TDI observablei = J2 d T^ m (15) 



m=l 



where 



d i = m max — uple polynomial of p variables Eij(t) (16) 

is the l-th TDI generator of a chosen generating set with I = l..l max . The variable l max 
is the total number of generators in the chosen generating set and m max is the number 
of data-flow variables considered in the vector e . The size, n with n < l ma x, of the 
smallest generating set of polynomials that cancel LF and OB noises, as well as m max 
or p depend on the level of modeling of the LISA mission. 



Indeed, one distinguishes between 1st, 1.5th and 2nd generation TDI. 
The First Generation TDI assumes constant (in time) and symmetric time delays 
(ty = tji). This leads top = 3 time-delay operators, a minimum of n = 4 TDI generators 
and m max = 6 data-flow variables (or 9 if LISA's internal motions are considered). The 
set of lst-generation-TDI observables forms the first module of syzygies over a ring of 
p = 3 variables [5J. For ideal, perfect, identical (stable, accurate, without shift nor 
drift) clocks beating the time t aboard the three stations, if the time delays are 
known exactly, within the lst-generation-TDI assumptions, the LF (and OB) noises are 
exactly cancelled in lst-generation-TDI observables formed from the data-flow variables 
recorded at each spacecraft. 

The 1.5th Generation TDI still assumes constant time delays, but they are not 
reciprocal anymore (tjj ^ tji). Hence, the number of time-delay operators is doubled 
p = 6, while n = 6 and m max = 9. The 1.5th-generation-TDI observables still form 
a module over a ring of p = 6 variables [8]. Hence, within the 1.5th-Generation-TDI 
assumptions, for ideal perfect identical clocks beating the time t and exact knowledge of 
tij, the LF and OB noises are still exactly cancelled in 1.5th-generation-TDI observables. 

The 2nd Generation TDI relaxes the constant time-delay assumption (t^ = tij(t)) 
with respect to the 1.5th Generation TDI. Consequently, the time-delay operators 
(p = 6) do not commute anymore. A new module still must be found. In reference 
[T3] , the authors propose 6 new, more complex, generalized Michelson- and Sagnac-type 
2nd-generation-TDI observables. For those, the order in which the delays are applied 
matters. Those do not remove exactly the LF and OB noises, even for ideal perfect 
identical clocks beating the time t and exact knowledge of t^, but those bring them 
down to an acceptable level for LISA. 
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Let us first consider the different TDI generations with respect to LISA geometry 
(orbitography and laser links) modeling. 

The first-generation-TDI assumptions are only met by a rigid, motion-less (this implies 
that no gravitational bodies are around to cause any motion) LISA constellation model. 
At this level of modeling, each time delay is given by the corresponding constant 
interferometric arm-length, slightly different from the generic arm-length L: = Ly/c 
with Lij = ~~x j — ~x V 

The 1.5-generation-TDI assumptions allow for a rotating (around its center of mass 
and around the barycenter) rigid LISA. At that stage, the Sagnac [9] and aberration 
effects cause the non reciprocity of time-delays pj. Indeed, for a same path-length, light 
rays travelling clockwise and counterclockwise do not take the same time (Sagnac effect). 
Additionally, there is a motion of the arm-length ij with respect to the barycenter 
causing an aberration effect. The model for LISA corresponding to 1.5-generation- 
TDI assumptions consists in classical motion (keplerian) of the three stations in the 
gravitational field of the Sun up to first order in eccentricity (eusA — 0.0096). 

The true LISA requires 2nd Generation TDI. A classical keplerian orbital motion 
around the Sun, if higher orders in eccentricity are considered, causes a flexing of the 
constellation: Lij(t) is a simple periodic function of t if only the Sun is considered, the 
so-called breathing of the triangle. If the newtonian perturbation of planets is consid- 
ered, the flexing becomes more complex. 

Moreover, if a gravitational relativistic description of photon time transfer is adopted, 
tij has additional time-dependant contributions due to gravitational relativistic effects. 
Indeed, in reference [21 Section III], a native, coherent, gravitational relativistic descrip- 

(0) (1/2) (1) 

tion of the laser link provides Uj = t ij + t ij + t ij, as a function of the positions 
and velocities of the emitting and receiving spacecraft at emission time. The 0th order 

in GM/c 2 ~ v 2 /c 2 , t ij, is the classical time taken by light to travel L{j(t) at velocity 
c, the l/2th order contains the Sagnac and aberration effects and the 1st order, light 

_ (0) (1/2) 

deflection or the so-called Shapiro delay. In reference [21 Section VI], these t ij, t 

(i) 

and t ij contributions to the photon flight time were evaluated numerically to 5 • 10 9 

m/c ~ 16.7 s plus a flexing of amplitude of 48000 km/c « 0.17 s, ~ 3- 10~ 3 s and < 10~ 7 

s respectively, assuming classical ephemerides of the three LISA stations. 

In reference |12| . relativistic ephemerides of the three LISA stations are provided. The 

authors compute numerically relativistic (t) — classical^) arid show that it reaches up 

_, (o) 

to about 3 km over a year, corresponding to an extra ~ 10 s in t ij. 

A realistic model of LISA's orbital motion also leads to geometric frequency shifts in 

expression (fl2"]l . In [2J Section V], the native, coherent, gravitational relativistic de- 

.. (l/2)« (3/2)*J 

script ion of the laser link also provides Aq = A G + A G + A g , as a function of the 
positions and velocities of the emitting and receiving spacecraft at emission time. The 
l/2th order is the Doppler shift due to the relative velocities of the emitting and receiv- 
ing stations in the special relativistic framework; it is reciprocal. Higher order terms 
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that contain gravitational relativistic effects such as the Einstein gravitational Doppler 

(i)« 

shift in Ac are non reciprocal. In reference [2J Section VI], these contributions to 
photon frequency shift were evaluated numerically assuming classical ephemerides of 

(1/2)^' (1)« (3/2)^' 

the three LISA stations, leading to A G ~ 7 • 1(T 8 , A G ~ 2 • 1(T 13 and A G ~ 2 • 10" 14 . 
Hence, the geometric frequency shift is irrelevant to LISA, owing to the mission's de- 
tection frequency bandwidth. 

Let us now consider the TDI analysis from the point of view of LISA time scales. 
TDI-observables (1st, 1.5th and future module for 2nd generation) should cancel exactly 
the LF and OB noises for ideal, perfect, identical clocks beating the time t aboard the 
three stations. However, general relativity taught us that the physical times scales 

12 3 

are the proper times scales r, r and r; that is the time beaten by the clocks aboard 
spacecraft 1, 2 and 3 respectively, if those are not constrained/synchronized^]. Strictly 

12 3 

speaking, r, r and r are different time scales which are not directly compatible, like 
apples and oranges. The natural common language between those time scales is a 
time coordinate t such as the one associated with the Barycentric Coordinate Reference 
System, called TCB. In the LISA mission, the events, which should thus be recorded in 
the corresponding proper time scale aboard a given LISA spacecraft k, are the different 
emission or reception event at that spacecraft. The locally measured frequency shifts 
(against the local proper frequency) should also be regularly recorded on board according 

to the local proper time scale (see Figured]). Each spacecraft k data-record should thus 

. k k 
contain T n em T n rec fc, 



r, at k with LF,OB,F,PM noises 



and 



A jfe (r) = A j7 



k 



t, at k with LF,Q,OB,PM noises + GW signal + G 



(17) 



(18) 



k 

where we adopted the notation r em / rec t for the n-th emission/reception event at space- 

k 

craft k on the time scale r. Consequently, to apply the TDI analysis, i.e. to compute 

u 

(jT5j) from (0 (TOl [U], Il2j) . we need to use equation (|7]) to convert from time scale r, in 
which the observed/simulated frequency shifts (TTT l [181) should be recorded, to time scale 
t. 



To summarize this section: a coherent general relativistic approach of TDI analysis 
is necessary. This means general relativistic modeling of laser links [2] (that is tij, since 
geometric frequency shifts can be omitted owing to LISA frequency bandwidth), LISA 
station orbitography [12] and time scales, which is the subject of the present paper. 

|| Indeed, equation ([7]) for stations k = 1..3 orbiting the Sun is similar to the GPS clock synchronization 
1 applied to clocks aboard the GPS constellation orbiting the Earth. 
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4. LISA Code: a LISA Simulator 

A LISA simulator aims at an intimate understanding of the mission, providing an 
overview of data analysis. This means that the simulator will be a precious tool to 
find and test strategies to detect GW, to test the specificities of interferometric config- 
urations for that purpose, to test the response of the detector (see for example LISA 
SIMULATOR [H, SYNTHETIC LISA p3] or LISA Code [llj). 

However, before getting to detection strategies, we have seen that the data must 
go through the TDI pre-processing analysis in order to remove (most of) LF and OB 
noises. A sine-qua-non condition for the LISA mission to work is for the TDI method 
to be efficient. 

Limitations on the effectiveness of the TDI technique come from the fact that the 
considered TDI generation is more or less adequate to the real orbitography and to the 
real laser links in the LISA mission. They also come from secondary noises (Q, F, PM, 
USO -Ultra-Stable Oscillator, that is the clock aboard LISA spacecraft-...) affecting the 
measurement, and certainly from the finite accuracy/precision of the quantities required 
to construct the corresponding TDI observables (see [EJ Section V]). 

Consequently, one aim of LISA Code is also to test TDI pre-processing efficiency 
for a realistic model of the LISA mission (orbitography, laser links, time scales) that is a 
coherent general relativistic rather than classical model, in presence realistic noises: LF, 
OB, Q, F, PM, USO. With that respect, the question that LISA Code addresses is: Is 
the residual LF plus OB noise left in a given TDI-generation observable acceptable? 
Presently, LISA Code uses a classical orbitography model (including constellation 
breathing plus rotation around its center of mass and around the barycenter), such 
as described in [6] and allows also for a static LISA or a rigidly rotating LISA. LISA 
Code adopts a general relativistic laser links description according to [2] to compute 
tij. LISA Code also contains LF, OB, PM and USO (shift and drift) noises; standard 
TDI observables of first and second TDI generation, as well as a procedure to compute 
arbitrary TDI observables. Each of the three USO noises initially represents intrinsic 
offset, difference in rate and drift of the clocks aboard one of the three LISA spacecraft. 
But via this USO-noise input in LISA Code, the proper time versus TCB time 

k 

transformations, A, described in Section [2721 can also be introduced. There will then 
be interesting questions to answer: How will the precision of LISA ephemerides, used in 

k 

LISA Code to compute both tij and A, affect the efficiency of a given TDI generation? 
And how will the constraints on LISA ephemerides accuracy derived with LISA Code 
compare with the theoretical constraints provided by Tinto et al. (about 16 m for arm- 
length [T4"l Section V A, (29) and (33)] and about 50 ns for clock synchronization [T5J 
Section V B, (39) and (40)]) for the 1.5th-generation standard TDI observables? 
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5. Conclusions 

We provided the relativistic analytical and numerical estimates for spacecraft 
proper times versus coordinate (BCRS) time transformation, for a classical orbitography 
in the gravitational field of the Sun. We showed that the difference between station 
proper time and BCRS time reaches about 0.5 s over a one-year mission and exhibits 
an oscillatory trend with a maximal amplitude of about 0.0014 s. The difference in rate 
of station proper time versus BCRS time is of the order of 1.5 ■ 10~ 8 . 
We then recalled how this time transformation fits in a general relativistic approach 
of Time Delay Interferometry, necessary to lower the laser frequency noise and optical 
bench noise down to the gravitational wave detection threshold. 
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Figure 1. The TDI method is developed in terms of the coordinate time scale, t, 
associated with the Barycentric Coordinate Reference System which is an artificial 
time scale (it can only be computed, not measured). Since each spacecraft i or j 
has its own physical time scale in terms of the proper time beaten by its clock r or 

r, a coordinate time t is the natural "common language" . However, events (such as 
emission or reception of a signal) and proper frequency shifts are recorded locally at 

i 

each spacecraft on their own proper time scale. Hence, a time transformation, A or 

3 . k k 

A, is required with A=t —t. 
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Figure 2. Differential proper versus TCB time transformation for spacecraft k (Eq. 
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Figure 3. Proper versus TCB time transformation (Eq. O for spacecraft k = 1, 2, 
3, integrated over a one-year mission. In this figure, for the sake of the illustration, 
there is an initial shift of 0.1; 0.2; 0.3 s for the clock aboard spacecraft k = 1, 2 or 3 
respectively. 
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time transformation (linear trend removed) (s) 




Spacecraft 1 

Spacecraft 2 

Spacecraft 3 

Figure 4. Differential proper versus TCB time transformation (Eq. [7]) for spacecraft 
k = 1, 2, 3, the linear trend is removed using a least-square fit method. 




Figure 5. Model for the LISA mission: double laser links between spacekraft i,j, k 
and double fiber links between the two lasers aboard each spacecraft. 
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